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We investigate the high frequency HaU effect on a two-dimensional triangular lattice with nearest- 
neighbor hopping and a local Hubbard interaction. The complete temperature and doping depen- 
dencies of the high-frequency Hall coefficient Rh are evaluated analytically and numerically for 
small, intermediate, and strong interactions using various approximation schemes. We find that Ru 
follows the semiclassical 1/qn* law near T = 0, but exhibits a striking T-linear behavior with an 
interaction- and doping-dependent slope at high temperature. We compare our results with previous 
theories as well as Hall measurements performed in the cobaltates. 
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I. INTRODUCTION 

The interpretation of the Hall resistivity is made dif- 
ficult even in relatively simple metals by the connection 
between the Hall constant and the relaxation time on 
the Fermi surfaccji In many cases the simple semiclassi- 
cal expression Ru = 1/qn* does not work and we must 
therefore deal with a transport measurement encoding 
much more information than just the effective density n* 
and the sign q of the charge carriers. 

In strongly correlated systems the Hall effect is even 
more difficult to interpret because interactions can have 
a large influence on the Hall resistivity. This influence 
moreover increases as the dimensionality of the systems 
decreases. There have been various attempts to describe 
the Hall effect in strongly correlated models with different 
geometries in 2D^ quasi-213^ and quasi-lD^i^ systems, 
but a general theoretical understanding is still lacking. 

Among the strongly correlated systems the triangu- 
lar lattice exhibits a unique property: it has the smallest 
possible closed loop with an odd number of steps (namely 
3). Anderson proposed that the model could have a spin- 
liquid ground state at commensurate fillings such as one 
electron per sitcjii These peculiarities make the triangu- 
lar lattice a very interesting system, which has been in- 
vestigated extensively in the last decades. In particu- 
lar, important differences between the Hall effect in the 
square and triangular lattices were pointed out in Ref. 0. 
Additionally, many studies have been motivated by the 
recent discovery of superconductivity in C0O2 layered 
compounds (cobaltates), which are good realizations of 
an isotropic 2D triangular lattice^ and in organic con- 
ductors of the BEDT family^ where one finds various 
structures resembling the anisotropic triangular system. 

Several Hall measurements have been undertaken in 
both organic superconductorsiSiii and cobaltates, espe- 
cially in the Naa;Co02 compound. In the latter, the 
anomalous linear increase of the dc Hall coefficienlii^ and 
a recent infrared Hall measurement^^ have motivated fur- 
ther theoretical work on this issu o^'^i^^'^^i^^ to investigate 



the effect of correlations in this system and their contri- 
bution to the Hall coefficient, but many questions re- 
main open. In the case of the organic conductors, the 
anisotropy present in the BEDT family complicates the 
problem further. 

In the present work, we study theoretically the Hall ef- 
fect in a 2D triangular lattice where electrons interact via 
an onsite Coulomb repulsion U. We calculate i?H in the 
high frequency limit^i^ where the probing frequency to is 
the largest energy scale of the problem, and we cover the 
whole range of interaction values using several approxi- 
mation schemes. 

The paper is organized as follows: in Sec. |TT]we intro- 
duce the model and the formalism used to compute the 
Hall coefficient Ru at high frequency on the triangular 
lattice. In Sec. IIIII we present our analytical and nu- 
merical results obtained in the whole (n, U, T) parameter 
space. We discuss in detail the doping and temperature 
dependence of i?H- Sec. [IV]is devoted to a discussion of 
our results and a comparison with other theoretical ap- 
proaches as well as experimental measurements. Finally 
our conclusions are given in Sec. |Vl and the subsequent 
appendices collect technical details. 



II. MODEL AND METHOD 

Our model is sketched in Fig. [TJ We consider an 
anisotropic triangular lattice with nearest-neighbor hop- 
ping amplitudes t and t' and an on-site Hubbard inter- 
action U . The Hamiltonian reads: 
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where cj, (cq) is the creation (annihilation) fermion op- 
erator, Ua is the fermionic number operator and {ij) are 
nearest-neighboring sites. The dispersion relation for this 
model (Fourier transform of iy ) is 

Ek = —2tcos{kxa) — it' cos{kxa/2)cos{kyaV3/2). (2) 
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FIG. 1: Top: Two-dimensional triangular lattice, a is the 
lattice parameter, t and t' are the hopping amplitudes for 
bonds along the x direction and for ±60° bonds, respectively. 
The current / flows along the x axis, the magnetic field B 
is applied along the z axis, and the Hall voltage is measured 
along the y axis. Bottom: Non-interacting density of states 
of the model in the cases t = t' — —1 and t = 2t' = 1. The 
DOS generically presents two van Hove singularities, which 
are degenerate in the isotropic lattice. The energy position 
of the van Hove singularity for t — t' — —1 corresponds to a 
band filling of 1/2 electron per site (1/4 filling). 



The correponding density of states (DOS) exhibits two 
van Hove singularities which are degenerate when t = t' 
(see Fig. [ij. UnUke in the square lattice, the DOS has 
no particule-hole symmetry, irrespective of the value of t 
and t' . 

We assume that a current / flows along the x axis 
and a dc magnetic field B is applied along z, hence 
a Hall voltage develops along the y axis (see Fig. [1]). 
In order to represent the magnetic field and the ap- 
plied ac electric field along x, we use the vector poten- 
tial A = A^'^^ -\- j4°', where for the magnetic part we 
choose the Landau gauge, A™^^ = Bxy, and A"^ de- 
scribes the electric field. The coupling between the lattice 
fermions and the electromagnetic field induces a Peierls 
phase in the hopping amplitudes which change according 
to tij -^ tij cxp(— ie /."' A ■ dl). 

The operator for the total current, J^ = j drj^{r), 
and the diamagnetic susceptibilities Xai(0) of th^ system 
are defined as usual: 
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where S = a?\/?>/2 is the unit-cell area, NS is the total 
system surface, (rifc) = {<^kCk) is the distribution func- 
tion and the thermodynamic average (■ ■ ■ ) is taken with 
respect to the Hamiltonian of Eq. ([T]). 

Performing the functional derivatives we find for the 
components of the currents 



J^ = eapt^4^Cfc^sin(fc^a) 
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where we have defined the vector r/ — [rj, 0) with 77 — 
\/3eBa/2. The diamagnetic susceptibilities resulting 
from Eqs (g]) and dH) are: 
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The Hall coefficient is defined as the ratio of the Hall 
resistivity to the applied magnetic field, Rn = py^/B, 
the Hall resistivity pyx being related to the conductivity 
tensor a^j^^ through 



Pyx 
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As shown in Appendix \^ it is possible to rewrite i?H as 
a high-frequency series where the infinite-frequency limit 
reads 
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and the remaining contributions are expressed in terms 
of a memory matrixi^iii Eq. ([8]) was originally derived in 
Ref.Uusing a different frequency expansion. Rui^ -^ 00) 
is expected to provide the dominant contribution at any 
finite frequency. The memory matrix formalism allows 
in principle to go beyond the infinite frequency approx- 
imation and compute corrections at finite frequency.— i^ 
It leads, in particular, to corrections due to interactions 
that vanish identically if C/ = 0. These corrections do 
not affect the sign of i?H which is entirely determined by 
i?H(w -^ 00). In the following we shall consider only the 
infinite- frequency contribution to Rn, Eq. ([5|), and adopt 
the notation R-ii{uj —> 00) = i?H- 



Strictly speaking, our results are valid provided the 
probing frequency is larger than any other energy scale 
in the system, to > inax{U,t,T}. The last two con- 
ditions, Lu > max{i,T}, are easily fulfilled experimen- 
tally in known triangular compounds, while the condition 
LU > U is more problematic. However, as we will discuss 
in Sec. IIVI in certain limits our results coincide with 
those obtained in Ref. [l3| under the opposite assumption 
cj <C [/, showing that this condition is not stringent. 

In order to evaluate Eq. ([5]), we calculate the commu- 
tator [Jx,Jy] from Eq. ([5]), and we use the diamagnctic 
susceptibilities of Eq. ^ to arrive at 
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As can be seen from Eq. ([9]) the high frequency Hall coef- 
ficient depends only on the distribution function (uk) as 
well as some geometrical factors. The interaction term 
in Eq. ([T]) therefore only influences i?H through its effect 
on (nfc). Since (uk) depends relatively weakly on U, we 
also expect the [/-dependence of Rn to be weak. Another 
implication of Eq. ^ is that at low temparature the be- 
havior of i?H can be interpreted in terms of an effective 
carrier concentration, as in the non-interacting case. 



III. RESULTS 

In the following we evaluate Rn in the whole domain 
of interaction values U with respect to the bandwidth 
W = 9t oi the system, by using four differents ap- 
proaches: exact calculation at [/ = 0, a pcrturbative 
expansion of the self-energy at [/ < W, a local approxi- 
mation to the self-energy, treated with dynamical mean 
field theory (DMFT) at U >W, and finally the atomic 
limit of the self-energy at U ^ W. 



A. Non-interacting case 

In the non-interacting case there are various limits in 
which we can obtain analytical results for i?^ {i.e i?H 
at [/ = 0): at zero temperature and band fillings near 
n — and n = 2, and at high temperature T :s> W. For 
intermediate fillings and temperatures, we compute R^ 
numerically by performing the sum in Eq. Q on a dense 
(2048 X 2048) discrete fc-point mesh. 



1. Zero temperature 

Here we restrict for simplicity to the isotropic case 
t' ~ t and we set the lattice parameter a ~ 1. Close 
to the band edges we can expand the various integrands 
of Eq. (jlOp and thus perform the k integrals. 

Near the bottom of the band the Fermi surface is made 
of two nearly circular electron pockets around (%^,0) 
and (3^, ■^)- In each pocket we have ^/c = Efc — yU ~ 
3t — |tfc^ — ^, where k is the momentum measured from 
the pocket center, and therefore fcp = |('^ ~ ^^/'t)■ The 
corresponding electron density is n = kp/ir. Writing sim- 
ilar expansions of A^, B^, and Ck close to the pocket 
center and performing the Brillouin zone integrations, we 
obtain the non-interacting Hall coefficient at low electron 
density: 



i?^(T = 0) = - 
ne 
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At sufficiently low density we recover, in the above ex- 
pression, the classical result R^ — 1/ne. 

Near the top of the band the Fermi surface is a nearly 
circular hole pocket centered at fc = (0,0). Close to 
this point we have ^k ~ ^6f + §tfc^ — /i, and therefore 
fcp = I (6 -I- fj,/t) . The corresponding density is obtained 
by substracting the contribution of the hole pocket from 
the maximum density: rih = 2 — fc|/27r. Similarly, for the 
functions Ak, Bk, and Ck we have to substract the con- 
tribution of the hole pocket from the contribution of the 
whole Brillouin zone, which turns out to be zero because 



^ Afc = ^ Bfe = ^ Cfc = 0. 
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Thus, for low hole densities n^ = 2 — n we find that the 
non-interacting Hall coefficient is given by 



KiT = 0) 



1 



1 



m 



-OK) 



(13) 



and as rih —> we have i?^ = —l/rihe. 

The complete density dependence of R^ calculated nu- 
merically at zero temperature from Eq. ^ is displayed 
in Fig [2] and compared to the limiting cases Eqs (fTTj) 
and (fT3|) . It is clear from this figure that the infinite- 
frequency i?H follows the well-known dependence of the 
dc Hall coefficient R-ai^u = 0) on the carrier charge den- 
sity. This indicates a weak frequency dependence of the 
non-interacting Hall coefficient at zero temperature, since 
the dc result is recovered from the infinite frequency limit 
of i?H- Furthermore this suggests, as we will discuss 
in more details below, that the frequency dependence 
should not be too crucial, even in the presence of inter- 
actions, for most band fillings. At U = the sign of the 
Hall coefficient is entirely given by the sign of the carri- 
ers, and it can be seen from Fig. [2] how the sign changes 
when the Fermi energy crosses the van Hove singularity 
of the DOS, and the Fermi surface shape evolves from 
electron to hole like. 
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FIG. 2: Non-interacting Hall coefficient R^ at zero tempera- 
ture as a function of the electron density n, for an isotropic 
triangular lattice with t' — t = ~1. The dashed lines indi- 
cates the clasiscal behavior at low electron and hole carrier 
densities. 



2. High temperature 
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FIG. 3: Temperature and electron density dependence of the 
non- interacting Hall coefficient i?^, for the isotropic trian- 
gular lattice with t = t' = —\. The dashed line shows the 
asympotic behavior described by Eq. (|16|l . Inset; Temper- 
ature and density dependence of the chemical potential /i, 
illustrating the relation /x ~ T at high temperature. 



If T 3> ^ the distribution function (rifc), which reduces 
to the Fermi distribution at [/ = 0, can be expanded 
in power of /3 = 1/T. This expansion must be done at 
constant density n, which requires that /3/z remains finite 
as /3 — > 0, in other words /.t ~ T at high temperature. 
Taking this into account we can deduce the relation be- 
tween /x and n, exp(— /3/i) = 2/n— 1, and write the Fermi 
distribution as 



K^k^ 
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Due to Eq. p2|) the fc-indepcndent terms in Eq. (|T4l) do 

:es the form: 



not contribute to i?^, which thus takes the form 
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Performing the BriUouin zone integrations we obtain 



K{T » t) = 
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This result is plotted in Fig. [3] together with the numeri- 
cally calculated full temperature and density dependence. 
The most striking feature of Eq. (fTB|) is the linear increase 



of i?j^ with T . The same linear behavior was obtained 
in Ref. |l8| at a; = 0, indicating a weak frequency depen- 
dence of i?j\ at high temperature. Our result shows that 
the T-lincar dependence of Ru is not due to interactions 
but to the peculiar topology of the triangular lattice. The 
sign of i?y at high T is determined by the sign of i, ir- 
respective of the density (see Fig. [3]). We attribute this 
property to the fact that at high enough temperatures the 
full band contributes to the Hall effect; hence the sign of 
i?H reficcts the dominant nature, electron or hole-like, of 



the band. As is clear from Fig. [H for i < the band is 
dominantly hole-like, while for t > it is electron-like. 

The relevance of result (fT6)) is that even without in- 
teractions in the system, the Hall coefficient has a lin- 
ear dependence at high temperature due to the geom- 
etry of the system, emphasizing the peculiarity of the 
triangular lattice. By contrast on the square lattice 
the same analysis yields a T-independent non-interacting 



7?o - 



n(2- 



at high temperature. 



B. Weakly interacting regime 

When interactions are present, the distribution func- 
tion (nfe) can be expressed in terms of the one-electron 
self-energy I](A;,ia;„) asji^ 
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with ujn = (2n+ l)7rT the odd Matsubara frequencies. In 
the weak coupling regime U < W, we evaluate the self- 
energy using conventional perturbation theory in U and 
we keep only the lowest order contributions of order U'^. 
For a local interaction like the Hubbard term in Eq. ([1]) 
there is only one diagram which is drawn in Fig. 31 The 
standard diagrammatic rules yield the following expres- 
sion for the self-energy: 



U' 



S(fe,zc^„) = -— ^ 



(18) 



/(Cfc2) [/(Cfci) - /(^fc+fci-fca)] ~ /(^fci)/(-Cfc+fci-fc2 
iiVn + $fci - S,k2 — ^fc+fci-fc2 

where f{S.k) is the Fermi distribution function. 



The numerical evaluation of Eq. p8l) is demanding due 
to the double momentum integration. This is particu- 
larly time consuming because our calculations are done 
at fixed density, and thus require to calculate E(fc, iuJn) 
many times in order to determine the chemical poten- 
tial. However it turns out that the momentum depen- 
dence of E(A;, iujn) in Eq. P^ is weak. This is illustrated 
in Fig. m where we plot the Brillouin zone average of 
the self-energy T,{iujn) as well as its standard deviation. 
The weak momentum dependence allows us to compute 
S(fc,za;„) on a coarse (typically 16 x 16) fc-point mesh, 
and then to interpolate using splines onto a dense mesh 
for the evaluation of (nfe) and eventually Rn. The Mat- 
subara sum in Eq. (fT7|) also requires special attention: 
the formal regularization of the divergence through the 
exponential factor is not suitable for a numerical evalua- 
tion of the sum. Wc therefore rewrite Eq. p7|) as 
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The LUn sum is now convergent and can be efficiently cal- 
culated via the truncation at some large frequency and 
the analytical evaluation of the remaining terms using an 
asymptotic expansion of the self-energy. 

The i?H resulting from perturbation theory are valid in 
the regime U < W <S^ui, with W = 9|t| the bandwidth of 
the system. As already anticipated the effect of a small 
U on the distribution (uk) is a subtle broadening, and 
as a result the dependence of Rn on U is very weak at 
low U. Fig. [5] provides an illustration of this weak de- 
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FIG. 4: Brillouin zone average of the real and imaginary parts 
of the self-energy Eq. (|18|) at low temperature T = 0.1, cal- 
culated using a 64 X 64 fe-point mesh (solid lines). The small 
standard deviations a (shaded curves) illustrate the weak mo- 
mentum dependence of E(fe, icOn). The dashed lines show the 
local self-energy resulting from the DMFT calculation (see 
Sec. line l|l . The density was set to n = 1.54, which is the 
value for Nao.7Co02 (see Sec. lIVp . Inset: Feynman diagram 
corresponding to Eq. (fT5)) . 



pcndence. As a consequence the non-interacting results 
of Sec. IIII Al are expected to give a fairly good account 
of the Hall effect for an interaction strength smaller than 
the bandwidth W. 

An important observation which wc can make from our 
pcrturbative calculations is that the momentum depen- 
dence of the self-energy is very small, i.e. the self-energy 
is almost local in real space. This suggests to approach 
the strong-coupling regime U >W hy assuming that the 
self-energy is exactly local. In the following section we 
study such local approximations to the self-energy, and 
we compare them to the result of the perturbation theory. 



C. Strongly interacting regime 

Assuming that the self-energy is local in first approx- 
imation, we investigate here two models for S(iw„) and 
their implications for the Hall coefficient i?H- The first 
approach is based on the dynamical mean field theory 
(DMFT)2S and requires to solve a difficult self-consistent 
quantum impurity problem. Due to numerical difficulties 
this method cannot be pushed to very high interactions 
and/or very low temperature. Our second approach is 
based on a simple analytical form for E(zcl'„), which is 
expected to be valid at U ^ W, and allows us to express 
(nfc) analytically in this limit. 



1. DMFT 

In the DMFT framework the local self-energy is ex- 
pressed as: 



Il{iuJn) == Qo'^i'i^n) - Q ^{i^n) 



(20) 



where Go is an effective propagator describing the time 
evolution of the fermions in the absence of interaction, 
and G is the full propagator, which takes into account 
the local Hubbard interaction. The calculation of Q from 
a given ^o amounts to solve the problem of a quantum 
impurity embedded in a bath. We do it by means of 
the quantum Monte Carlo Hirsh-Fye algorithm^! as de- 
scribed in Ref. |20| . From the requirement that G coincides 
with the local Green's function of the lattice, i.e. 



G{iio„ 
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(21) 



one can deduce the self-consistency condition 

Go^iiuJn) = l/D [iLUn - T.{iuJn)] + T,{iuJn), (22) 

where D{z) = / d^ D{^)/{z — ^) is the Hilbert transform 
of the DOS D(^) corresponding to the triangular lattice 
and shown in Fig. [TJ Once the self-consistent Goi'^'^n) 
is obtained, the corresponding self-energy T,{iujn) is in- 
jected in Eq. P^ to compute i?H- 
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FIG. 5: Evolution of the high-frequency Hall coefficient with 
U calculated, for an isotropic triangular lattice t — t' — ~1, 
using different approximations at T = \t\ and ?i = 1.54 



In Fig. [4] we compare the DMFT self-energy with the 
Brillouin zone average of the pcrturbative expression 
Eq. (|18|1 . both calculated at t/ = 4. It can be seen that 
the frequency dependence and the order of magnitude of 
the two quantities are very similar, suggesting that the 
self-energy is dominated by the U^ term and therefore 
the domain of validity of the perturbation theory is not 
limited to very small U. On the other hand it shows that 
the DMFT, although it is not a pcrturbative approach, 
provides a smooth transition from the weak to the strong- 
coupling regimes. This is further illustrated in Fig. [5] 
where we see that the values of i?H calculated by pertur- 
bation theory and DMFT coincide up to C/ « 4|t|. At not 
too low temperature the DMFT calculation is reliable up 
to interaction strengths comparable to the bandwidth W. 
We have performed DMFT calculations at U > W, but 
since these results could be affected by systematic sta- 
tistical errors in the Monte-Carlo summation, they are 
not shown in Fig. [5l At t/ ^ VF it is expected that the 
DMFT result approaches the atomic limit in which ac- 
curate calculations can be performed, as discussed in the 
next paragraph. 



2. Atomic limit 

In the limit of very strong interactions U ^ W we 
assume that the self-energy approaches its atomic limit 
given by the expression (see Appendix [C]) : 






n/2(l-n/2)[/2 
, + ^at - (1 - n/2)U 



(23) 



with /Xat the chemical potential in the atomic limit, not 
to be confused with the lattice chemical potential fj.. Us- 
ing this expression in Eq. ([T7)) it is possible to evalu- 
ate analytically the sum over Matsubara frequencies and 
thus to obtain a closed expression for (n^) (Appendix 



[C|) . In Fig. [5] we show the Hall coefficient calculated with 
the atomic limit of the self-energy in the whole range of 
interaction values. Ru obviously converges to the non- 
interacting limit at low U since the atomic self-energy 
vanishes at C/ = 0, and provides a good interpolation 
between the weak and the strong-coupling regimes. At 
intermediate values U '^ W the atomic limit is not reli- 
able, although it gives the correct order of magnitude for 
Ra . Fig. [5] also shows that i?H saturates at sufficiently 
large U. 

In Fig. [6] we display the temperature and density de- 
pendence of i?H at U ~ 500|<| , which is a value typical for 
the cobaltate compounds as discussed in the next section. 
We have selected four densities corresponding to the bot- 
tom and top of the lower and upper Hubbard bands (see 
also Fig. [S] below). Like for U = wc find a T-lincar 
increase of i?H at T > W. Due to the Mott gap, how- 
ever, the density dependence of the slope is not the same 
as for U = 0. The slope can be obtained explicitly by 
sending U to +oo and performing the high-temperature 
expansion as in Scc. lIII A2l (see Appendix [C|) . The result 
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very similar to Eq. p6p except that the slope ex [4(5(1 — 
S)]'-^ replaces [27i(2 — n)]^^, where S ~ |'^ ^ 1| measures 
the departure from half-filling. The U = oo result of 
Eq. ([M)) is displayed in Fig. [51 and correctly describes our 
high-temperature results at U ~ 500|t|. The differences 
observed at n = 1.05 in Fig.[n]reficct the fact that close to 
half-filling the slope of the high temperature Eh depends 
strongly on the interaction and is not saturated even at 
U = 500|i| (see also Fig. 7). Away from half-fiUing the U 
dependence of the slope is weaker, and Eq. (24) is valid 
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FIG. 6: Temperature and density dependence of the high 
frequency Hall coefficient calculated in the atomic limit at 
U = 500|f|. The dashed line shows the asymptotic behavior, 
Eq. p4|) . and 5 = \n — 1\. For densities close to half-filling 
(n = 1.05), Rn(T) deviates from the asymptotic behavior (see 
text). 



for lower interaction values. 



IV. DISCUSSION 

The various approximations presented above allow us 
to calculate the Hall coefficient on the triangular lattice 
for all interactions strengths U and all temperatures T. 
The main limitation of our approach, in view of a com- 
parison with experimental systems, is that our results 
are in principle valid in the limit W^U <^ lo, because 
they are based on a high-frequency expansion. The first 
criterion, W -^ ui, is not too difficult to satisfy for re- 
alistic compounds if the measurement of the Hall effect 
is performed at optical frequencies. The second crite- 
rion, U <S^ u, seems more problematic since interaction 
strengths can be as large as several eV, at the upper 
edge of the mid-ultraviolet frequency domain. However, 
we have seen (Fig. O that at [/ = and T = the Hall 
coefficient calculated at lo ~ oo coincides with the uj ~ 
dc value, and at [/ = and T ^ t,we obtained the oj = 
results of Ref. [iS- AH this suggests that the frequency 
dependence of i?H is weak in the non-interacting case. 

At the other extreme of the parameter space, U = oo 
and T ^ W, we can compare the result of the atomic 
limit approximation with the result of the t-J model J^ In 
the latter model U is considered infinite from the outset, 
so that the high-frequency and high temperature expan- 
sion of Ref. [ij is in fact valid at frequencies uj < U. We 
plot in Fig. [7| the density dependence of i?H obtained in 
both models at [/ = c» and T > W. The small quanti- 
tative difference between the atomic limit at [/ = oo and 
the t-J model shows that these two ways of treating the 
U = oo limit are not equivalent: they differ, in particular, 
in the renormalization of the kinetic energy by the inter- 
action. However the two models give, where they can 
be compared, very similar behaviors. This reinforces the 
idea that the frequency dependence of i?H is weak. Ex- 
act diagonalization on small clusters also indicate such a 
weak frequency dependences^ This strongly suggest that 
our results could also be valid at w < C/, and therefore 
be relevant to interpret experiments performed in this 
regime. The atomic-limit approach has the advantage to 
give access to the full temperature dependence (Fig. [5]) 
as well as the [/-dependence as shown in Fig. [71 while the 
calculation of Ref. [Ij is valid at t/ = oo and T ^ W. 

The evolution of i?H with temperature is of particu- 
lar interest since it is most easily probed experimentally. 
A linear increase of Rn with temperature, without sat- 
uration at high T, was reported in Ref. |l3| for the t-J 
model. Our results show that the Coulomb interaction 
is not responsible for this effect which is also present at 
U = (Fig. [3]) and is therefore a consequence of the pe- 
culiar geometry of the triangular lattice. However the 
interaction controls the density dependence of the slope 
which changes smoothly from [2n(2 — n)]~^ at [/ = to 
[4(5(1 — 5)]~^ at U — oo. This is further corroborated in 

Fig.m 
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FIG. 7: Density dependence of the Hall coefficient within the 
T-linear regime ai U — and in the atomic limit at U ^ |t| 
(solid lines), compared to the result of the t-J model (Ref. Il3l. 
dashed line). 



The sign of i?H turns out to be independent of n and 
U at high temperature, unlike in the square lattice where 
i?H changes sign at n ~ f . The situation is different at 
T = 0. In the non-interacting case Rn changes sign at 
quarter filling and can be simply interpreted in terms of 
the carrier density (Fig. ^. We have also investigated 
the T = density dependence of i?H at large U as shown 
in Fig. [8l The interpretation in terms of the carrier dcn- 
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FIG. 8: (a) Density dependence of the Hall coefficient at 
T = Q and U — 500 |f| calculated in the atomic limit approx- 
imation, (b) Density of states for various electron densities, 
showing the lower (LHB) and upper (UHB) Hubbard bands. 
The shaded regions indicate the occupied states and the po- 
sition of the chemical potential. 



sity remains qualitatively valid, provided one takes into 
account the splitting of the DOS into the lower and up- 
per Hubbard bands. These two bands are displayed in 
Fig. [51d, where it can also be seen that the DOS keeps 
qualitatively the same shape as for U = 0, but the width 
of each band varies strongly with the density n. Due 
to this band renormalization the sign change of Rn at 
n < 1 does not occur at quarter filling, but a little below. 
Comparing Fig. [8^ with Fig. [7] one easily understands 
why the temperature dependence of Rn has to be more 
pronounced slightly above n = and n = 1 than slightly 
below n = I and rt = 2, as can be also seen in Fig. [S] 

Among the Hall measurements reported for layered 
compounds with a triangular lattice structure, there is 
one performed at finite frequency by Choi et alM- on the 
cobaltate Nao.7Co02. This material is composed of two- 
dimensional layers of edge-sharing CoOe octahedra sepa- 
rated by an insulating layer of Na"*" ions, leading to a tri- 
angular lattice of C0O2 units32, ARPES measurementsSl 
indicate that the triangular lattice is isotropic with an 
estimated hopping amplitude of < = —10 meV and an 
effective Hubbard energy U ^ b eV. From the radius 
of the Fermi-surface hole pocket observed in ARPES, 
/cp = 0.65 ± 0.1 A"-'^, we deduce an electron density 
n = 1.54. Choi et al. measured the temperature de- 
pendence of both the dc and ac Hall coefficients up to 
room temperature. The ac measurement was performed 
at a; = 1100 cm^^ « 12|i|- The experimental conditions 
thus satisfy T,W <uj <^U. 

We note, however, that there are discrepancies between 
different sets of experimental dataJ^i^ The behavior of 
the dc i?H above T = 250 K is consistent with the lin- 
ear increase predicted by the various theoretical models. 
By adjusting this models on the dc experimental data 
at high temperature (dotted line on Fig. ^ we obtain 
an independent determination of the hopping amplitude 
t, namely t = —7.4 meV using the atomic limit model 
Eq. ([241) and t = -5.7 meV using the t-J model. This 
values are in good agreement with the ARPES determi- 
nation of t. 

The organic conductors of the BEDT family present 
several compounds with an anisotropic triangular struc- 
ture. Unfortunately we are not aware of any measure- 
ments of the ac Hall effect which we could compare to 
our calculations, although measurements have been done 
at zero frequency in these materials J^iii 
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FIG. 9: Comparison of the Hall coefficient in Nao.7Co02 as 
measured by Choi et ali^ at oj = (empty circles) and uj = 
1100 cm^^ (full circles) with available theoretical models. The 
blue solid line is the high-frequency result in the atomic limit 
and the dashed line shows the high-temperature result in the 
t-J model.— Theoretical curves are calculated aX t — t' = 
— 10 meV, U — 5 eV, and n = 1.54. The dotted line shows 
the best fit of the experimental data at high temperature with 
theoretical models (see text). 



These conclusions apply provided the probing frequency 
is larger that the other energy scales of the problem, and 
that the electron self-energy remains essentially local for 
strong interactions. 

Although we do not expect that a possible momentum 
dependence of the self-energy can have a strong effect on 
i?H, and we have argued that the frequency dependence 
of i?H is probably weak, it is clear that for understand- 
ing the anomalously large i?H(^) measured experimen- 
tally in Nao.7Co02 in the mid-infrared range, one would 
have to extend the approach in order to cover the do- 
main of intermediate frequencies. Another possibility is 
that the simple one-band model considered in this study 
would not suffice to capture the detailed properties of the 
materialsi^ Experiments conducted as a function of w, 
as well as measurcinents of other materials with a trian- 
gular structure, would be very helpful to elucidate the 
peculiarities of the Hall effect in triangular compounds. 



CONCLUSION 



The theoretical high-frequency Hall coefficient in the 
two-dimensional triangular lattice exhibits two different 
characteristic behaviors at low and high temperatures: 
near T = 0, i?H resembles the classical dc Hall coefficient 
l/qn* where q and n* are the carrier charge and density, 
respectively; at temperatures higher than the bandwidth, 
on the other hand, i?H shows a remarkable T-linear be- 
havior with a density- and interaction-independent slope. 
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APPENDIX A: MEMORY MATRIX FORMALISM 

The use of the memory matrix aUows one to perform 
finite-order frequency expansions, which are singular for 
the conductivities due to their resonance structurCfii 
This approach has been used in previous works to study 
transport properties in Luttingcr hquids;^ as well as the 
Hall effect in the 2D Hubbard modcl^ and in quasi one- 
dimensional systems i^ 

As we want to calculate the Hall resistivity pyx, we 
start from the general relation between py^ and the con- 
ductivity tensor u^^: 



Pyx — 



'xy 



(Al) 



Then we rewrite the conductivity tensor in terms of a 
memory matrix M{uj) as^^ 

cT^iu;) =i{u;t + x(0) [^ + iM[u)] X-'(O)}"' x(0) 

(A2) 
where a^ denotes the transpose of cr, Jl is called the fre- 
quency matrix, and M{u}) is the memory matrix. x(0) 
is a diagonal matrix composed of the diamagnetic sus- 
ceptibilities in each direction, XtJ.v{^) = <^/ji/X/i(0). The 
frequency matrix Jl in Eq. (|A2p is defined in terms of the 
equal-time current-current correlators as^ 



^ ^tiu 



sx,m 



{[j.^jA) 



(A3) 



with S the sample surface (in a two-dimensional system) . 
Now we invert Eq. (|A2p and express the memory matrix 
M in terms of the conductivity tensor. For the Hall 
coefficient i?H we need only the off-diagonal term Mxy 
given by 



iMxy{uj) 



iXy{0)<^xy{uj) 



(7xx{i^)cryyiuj) + alJuj) 



n. 



(A4) 



This implies that the Hall coefficient i?H = Pyx/ B can be 
expressed in term of the frequency and memory matrices 
as 



Rai'^) = 77TT 1™ „ 

IXyiQ) B-.0 B 



(A5) 



Since the memory matrix vanishes as uj^"^ at high fre- 
quency, we see that the Hall coefficient is given by Eq. ([5]) 
in the infinite frequency limit. 



time slice in order to single out the Hubbard interaction. 
In a second step the interaction is decoupled via the in- 
troduction of an Ising variable in every time slice. The 
Green's function Q{t) is finally calculated by averaging 
over the ensemble of configurations of the Ising variables 
using a Monte-Carlo sampling and local updates. In our 
calculations at /3 = 1 and [/ ^ 20 we take L = 128 and 
we keep 10^ out of the ~ 10® configurations visited. The 
numerical accuracy of the calculated Q{t) is estimated to 
be ^ 10""^ at the highest U values, and closer to 10"''^ at 

In order to calculate the self-energy and solve the 
DMFT self-consistency condition, Eq. ([22)) . we need to 
Fourier transform Q{t) from imaginary time to imaginary 
Matsubara frequencies iw„. In traditional implementa- 
tions of the algorithm this step is performed through a 
cubic spline interpolation of G{t). Because cubic splines 
arc non-analytic, however, the resulting Fourier series are 
unreliable at frequencies above ~ L/ (3. Instead of an in- 
terpolation, we have performed a fit of Q{t). Our fitting 
function is a discrete form of the spectral representation, 
Q{t) = — J de A{e)e^°'^ f{—e), which we express as 



M 



g{T)^-Y,A,e-'^-I{-e,) 



(Bl) 



M 



1 Aj 



1. The number AI of poles 
are determined by adding more 



with Aj > and Y. 
£j and their weight Aj 
and more terms in Eq. (jBip . until the fitted function 
matches all QMC data points within a numerical toler- 
ance, which we take as the estimated accuracy of G{t). 
The Fourier transform is then simply 



Q{i^n) 



M 

E 



A, 



lUJn 



(B2) 



The calculated self-consistent propagators G{t) and 
tJo(T) are displayed in Fig. [TOl together with the fits to 
Eq. (IBTI) . 

Solving Eq. (P^ at fixed electron density requires to 
determine the chemical potential p self-consistcntly. In 
our calculations we perform the search for both the self- 
consistent Qq and p in one shot using a global minimiza- 
tion procedure. As a result the self-consistent solution 
can be reached in typically less that 20 iterations. 



APPENDIX C: SELF-ENERGY AND 

DISTRIBUTION FUNCTION IN THE ATOMIC 

LIMIT 



APPENDIX B: 



CALCULATION OF THE DMFT 
SELF-ENERGY 



We evaluate the local self-energy in the DMFT frame- 
work using the Hirsh-Fye algorithm^i as described in 
Ref. [20. In this method the imaginary-time axis [0, j3[ is 
cut into L slices, and the Trotter formula is used in each 



In the case U ^ t one can treat the Hamiltonian 
Eq. ^ using a perturbative expansion in Uj/U. The 
atomic limit is the zeroth-order term of this development, 
and it corresponds to a collection of disconnected sites 
with four possible states on each site. This limit is not 
very useful since there is no hopping and thus no possi- 
ble transition below the Hubbard energy U. In order to 
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FIG. 10: DMFT imaginary-time propagators calculated at 
/3 = 1 and n = 1.54 for various interaction strengths. The 
symbols show the QMC results on the discrete-time mesh. 
The solid lines are the fit, Eq. ()B1|I . used to evaluate the 
self-energy shown in Fig. ID 



retain the low-energy dynamics of the problem we adopt 
an hybrid approach, where the free dispersion is used in 
the lattice Green's function together with the self-energy 
evaluated in the atomic limit. The atomic self-energy is 
obtained by diagonalizing the Hamiltonian Eq. ([T]) with 



0, which leads to the atomic Green's function 



Ga 



l-n/2 

iuJn + Mat 



n/2 



Mat 



L/' 



(CI) 



while the non-interacting Go. at — l/(*^n + Mat) results by 
putting U = Q. From Dyson's equation, S = Gq"^ — G~^, 
we deduce the atomic self-energy displayed in Eq. (^5)) . 
Here Mat is the chemical potential in the true atomic 
limit — i.e. the limit where the lattice Green's function 
takes the form (jCip . and therefore the electron density 
is given by n = (2 - n)/(-/Xat) + nf{U - Mat) with / the 
Fermi function. We can invert this relation and express 
/iat explicitly in terms of the electron density as 



Mat 



log 



1 



+ y(^-l)' + ^-^^(^-l) 



Using the atomic self-energy Eq. ()23p as an approxima- 
tion to the exact self-energy in Eq. (fT7|) , wc evaluate ana- 
lytically the lattice distribution function (n^). Let's first 
remark that 



1 



«W„ — ^fc — Sat(iWn) iUn - E^ iuj„ - E^ 



with 



E^ = (Cfc ± Afe + [/ - /iat)/2 

. ■^fc + Afc + [/ - /iat 
-^k = 1 X 

4Afc 



At 



(Cfc + Afc - £/ + Mat) (Mat -U + nU/2) + n^l^tU 

Cfc(Mat -U + nU/2) + nn^tU/2 
Vi^k + U + Mat)2 + 2(n - 2)(a + Mat)C/. 



As a result the Matsubara sum in Eq. pT|) is easily per- 
formed to yield 



{nk).t = Ak.fiE+) + il-Ak)fiE-^ 



(C2) 



Within this approximation it is also straightforward to 
perform the infinite U limit. Taking into account that 
both /i and /iat are either of order t (if n < 1) or of 
order U {ii n > 1) we find that Ak approaches n/2 as 
U increases toward +oo. Likewise, if n < 1 we have 



E7 



U and E: 



t while if n > 1 we have E"^ 



t and 



— [/. Hence we find 



/,T \U=cc 
("fc)at 



vl-f)/[(l-f)?fc-fMat] n<l 



r/[fa-(i 



Mat 



1-f n> 1 



where we have introduced /Iat = Mat ^ ^ ^ ^Ti log[(l ~ 
n/2)/ {n — 1)] and S,k ^ £k ^ fj- with fl = ^i — U . For the 
purpose of evaluating the high-temperature behavior of 
the Hall coefficient at t/ = oo, wc finally expand the dis- 
tribution function in powers of P following the procedure 
described in Sec. IHIA2I 



/^ \U=oo 



- 7i(l - n)e4 + 0(/?2) 



n< 1 



f -(n-2)(l-n)£4 + 0(/32) n>l 



Comparing with Eq. p4)) . which is valid at C/ = 0, we 
see that the only difference between the high-temperature 
behaviors of Rn at C/ = and C/ = oo is the n-dependent 
slope, and we easily deduce that 
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T/t 
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By introducing S = \n — 1\ which measures the doping 
with respect to half-filling, these two cases can be recast 
in one single expression shown in Eq. 
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